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The optical conductivity spectra of the rattling phonons in the clathrate BagGai6Ge3o are investigated in 
detail by use of the terahertz time-domain spectroscopy. The experiment has revealed that the lowest-lying 
vibrational mode of a Ba(2) 2+ ion consists of a sharp Lorentzian peak at 1.2 THz superimposed on a broad 
tail weighted in the lower frequency regime around 1.0 THz. With decreasing temperature, an unexpected 
linewidth broadening of the phonon peak is observed, together with monotonic softening of the phonon peak 
and the enhancement of the tail structure. These observed anomalies are discussed in terms of impurity 
scattering effects on the hybridized phonon system of rattling and acoustic phonons. 
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1. Introduction 

A decade or more ago, certain compounds of the 
clathrates family 1-3 were found to exhibit metallic elec- 
trical conductivity but a heavily suppressed thermal con- 
ductivity resembling glassy materials, which lead to their 
potential application for thermoelectric devices. In gen- 
eral, they have a network structure of atomic cages in- 
cluding guest atoms. These guest atoms inside the cages 
give rise to THz-frequency optical modes, which have 
been observed from the specific heat, 3-5 inelastic neutron 
scattering (INS), 6,7 and optical spectroscopy 8-11 mea- 
surements. Such a low-frequency mode is considered not 
only to interfere but also to hybridize strongly with the 
heat-carrying acoustic phonons of the cage network, and 
thus is expected to play key roles in various anomalous 
phenomena mentioned above. On the basis of such per- 
spective, extensive studies have been conducted on hy- 
bridized phonon systems and their related properties. 

The type-I clathrate BasGai6Ge3o, hereafter abbrevi- 
ated as BGG, has a cubic structure with space group 
Pm3n, presenting 2 dodecahcdral and 6 tetrakaidecahe- 
dral cages formed by a network of Ga and Ge. It is a 
so-called Zintl compound that, for the charge compen- 
sation, each cage encapsulates a guest Ba 2+ ion, labeled 
Ba(l) and Ba(2), respectively. In general, a fine tuning 
of the Ga/Ge concentration slightly off the stoichiome- 
try turns the system into a heavily doped semiconductor, 
and hence controls both the charge carriers sign and den- 
sity. 3 ' 12-15 

Each Ba 2+ ion is loosely bounded to the local poten- 
tial of an electronegative cage, the tetrakaidecahedron in 
particular, and therefore vibrates with larger amplitude. 
Such a localized vibrational mode of Ba(2) in the over- 
sized cage is called as rattling phonon and shows strong 
anharmonicity. This anharmonicity results in a so-called 



softening whereby the vibrational frequency decreases to- 
ward low temperature, as demonstrated experimentally 
by Raman scattering 9 and theoretically by a quasihar- 
monic approximation to a quartic term of the guest ions 
displacement . 1 7 

In addition to this anharmonicity, the hybridization 
with the acoustic phonon is a significant characteristic 
of the Ba(2) 2+ rattling vibrations. In 2008, INS mea- 
surements on BGG have clarified the phonon dispersions 
at room temperature. 7 A clear evidence is obtained for 
an avoided crossing between the rattling vibration of 
Ba(2) and the cage acoustic modes. In comparison to the 
zone-boundary acoustic phonons of ~ 2 THz, the rattling 
phonon with frequency of ~ 1.2 THz is indeed low enough 
to become hybridized in an avoiding manner such that 
both phonon modes are strongly coupled and no longer 
independent of each other. It was pointed out that such a 
hybridization can suppress thermal conductivities at low 
temperatures 7 ' 16 

Under thus background of research, we have been 
highly motivated to study the charge dynamics of these 
rattling vibrations, and so far we have reported the opti- 
cal conductivity spectra of BGG 10 and the isostructural 
BGS (Ba 8 Gai 6 Sn 3 o). 1:L For the vibration of Ba(2) 2+ ions, 
there- are two infrared-active normal modes of T± u sym- 
metry, labeled as v\ and perpendicular and parallel 
to the local fourfold inversion axis, respectively. 9, 10 For 
BGG, the spectral sharpening of the v\ mode associated 
with softening toward low temperature was reported in 
the previous work. 10 This result was ascribed to both 
unequally spaced vibrational levels in a single-well an- 
harmonic potential, or the on-center potential, and the 
Boltzmann factor, expected from the local anharmonic 
potential model. 18 For BGS, with decreasing tempera- 
ture, the single broad peak of v\ mode became split 
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into two subpeaks and also showed a linewidth broad- 
ening. The former result was consistently understood by 
assuming a multi-well anharmonic potential, or the off- 
center potential, 18 revealed by structural analyses. 5, 15 
However, the latter finding was quite unexpected be- 
cause the low-temperature spectra becomes sharp due 
to reduction of thermal effects in decay process in gen- 
eral and suppression of thermal excitation in the local 
anharmonic potential model. 18 So far this anomaly has 
had yet to be explained, although we have suggested the 
relevance of interactions with other excitations such as 
acoustic phonons and/or charge carriers. 11 

As mentioned above, the temperature dependence of 
the linewidth of the rattling phonon spectra seems to 
show opposite behavior between BGG and BGS; the for- 
mer shows sharpening, while the latter shows broadening 
with decreasing temperature. Then there arises a ques- 
tion whether those contrasting behaviors can be ascribed 
to the difference of the charge carrier and the potential 
shapes, i. e., on-center or off-center. It is noted that, from 
the spectra of p-type BGG, we can clarify the dynamical 
properties of the rattling phonon without taking into ac- 
count any effects of charge carriers, that is, we have the 
information on roles of anharmonicity in the linewidth. 
In fact the dc resistivity increases sharply by a few order 
of magnitude at low temperature, 14 indicating that the 
carrier density and hence the interaction with rattling 
phonons might be negligibly small. 

Looking back to previous BGG spectra (see Fig. 4 in 
Ref. 10 ), one finds some spectral weight appearing in the 
low energy region of the V\ peak, which may suggest 
a double-peak structure and may lead to a modifica- 
tion of the analysis. In order to have a definite answer 
on this question, the frequency resolution must be im- 
proved much more than in our previous works, where 
the linewidth and the experimental resolution were in 
the same order. Thus we study in more detail the optical 
conductivity of BGG by improving the resolution in fre- 
quency and by analyzing carefully a spurious interference 
effect of multiply reflected light within a sample. 

In the present paper we report newly obtained opti- 
cal conductivity spectra in the THz range, with an im- 
proved frequency resolution over the previous work. 10 
The data show that the lowest-lying spectra of the rat- 
tling phonons consist of a sharp Lorentzian peak at 
1.2 THz superimposed clearly on a broad tail around 
1.0 THz. This paper consists as follows. In §. 2, we 
describe the experimental procedure and demonstrate 
how the resolution is improved by using a thicker sub- 
strate. In §. 3, we show the optical conductivity spec- 
tra where the interference effect is quantitatively taken 
into account. From the temperature dependence of the 
peak frequencies and linewidths, we conclude that the 
linewidth broadening appears toward low temperature. 
In §. 4, we discuss these results by use of the strongly 
hybridized model of the rattling and acoustic phonons 
with impurity-type scatterings. After commenting that 
the anharmonicity works as the mean field to induce 
a frequency softening, we suggest that the strong low- 
frequency hybridization between rattling and acoustic 



phonons, softening of the rattling phonon, and impurity- 
type scattering of phonons can be causes of the anoma- 
lous linewidth broadening toward low temperature. Fi- 
nally, in §. 5 we summarize the present works. 

2. Experimental procedure 

Single crystals with p-type carriers were grown by a 
self-flux method. 4, 14 To obtain optimum transmitting 
signals, disks of about 5 mm in diameter were attached 
with adhesive on a sapphire substrate and polished down 
to d =10-20 /im in thickness. Terahertz time-domain 
spectroscopy (THz-TDS) measurements covering the fre- 
quency range of 0.2-3.0 THz (6.7-100 cm -1 ) were carried 
out with a Tochigi Nikon RT-20000 spectrometer, which 
uses a standard technique for the transmission config- 
uration. 10, n ' 19 The optical conductivity is numerically 
determined by using the equation taking into account 
multiple reflections in the sample and adhesive. Four 
samples in total were used both for determining the op- 
timum thickness that gives highest signal-to-noise ratio 
and for more carefully mapping an interference effect. By 
using a thicker sapphire substrate of 2.0 mm rather than 
0.5 mm in our previous work, 10 the frequency resolu- 
tion in the present work has been greatly improved from 
about 100 GHz to about 25 GHz. For further details on 
our time-domain spectroscopy, refer also to Refs.. 10, u ' 19 
See also Appendix A and B for detailed discussions of an 
interference effect and a frequency resolution. 

3. Experimental results 

Figure 1 (a) and (b) show the real part eri(aj) and 
the imaginary part 02 (w) of the complex conductivity, 
respectively, at temperatures 296 K, 150 K and 5.7 K. 
The cti(w) spectrum presents distinct peaks superim- 
posed on almost constant conductivity due to charge car- 
riers, whereas 02 (w) shows the corresponding frequency 
derivatives, as previously observed. 10 The low-lying spec- 
tra observed around 1.2 and 1.8 THz, assigned as v\ and 
V2, show a drastic change with decreasing temperature, 
associated with the phonon softening. In contrast, the 
higher frequency peaks ^3 (a collective cage mode) and 

(a Ba(l) mode) around 2.4 THz arc weakly dependent 
on temperature. 

As indicated by arrows in Fig. 1 (a), we have identi- 
fied a few odd structures most clearly seen in the 5.7 K 
spectrum; for example a hump around 1.5 THz and dip 
around 2.2 THz. Those small structures are attributed to 
an interference effect between multiply reflected lights in 
the present highly dielectric material (See also Ref. 10 ). 
The dotted line in Fig. 1 (a) shows the interference pat- 
tern, which is estimated as described in detail in Ap- 
pendix B. In the lower panel of Fig. 1 (a), the interference 
pattern is also shown by an expanded scale. Those spuri- 
ous effects are observed weakly but clearly for frequency 
regimes in between phonon peaks or at low temperatures 
where any carrier contributions are hardly traced. We 
should mention here that a broad sideband at the lower 
frequency side of the z/i peak still appears even after the 
background correction is made. Since such a structure 
is reproducibly observed in our measurements for other 
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Fig. 1. (color online) (a) Real parts cri{uS) and (b) imaginary 
parts (72 (co) of complex conductivity spectra of p-type BGG. Dot- 
ted line and lower panel in (a) display the estimated interference 
pattern as described in Appendix B. (c) and (d): Phonon spec- 
tra for rattling mode ui at 5.7 K and 296 K. Experimental data are 
shown by open symbols and the fitting curves with two Lorcntzians 
by solid lines, see the text. For convenience, 1 THz = 33.3 cm -1 = 
48 K = 4.14meV. 

samples with different thicknesses, we conclude that this 
remaining sideband structure is intrinsic. For detailed 
discussions, see also Appendix B. 

From the corrected spectra, the temperature- 
dependent, non-Drudc-likc contributions of charge 
carriers are subtracted by assuming a linear dispersion 
with some gradient as an adjustable parameter for the 
sum rule to hold for oscillator strengths. In Fig. 1 (c) 
and (d) we show the peak profiles for the rattling v\ 
mode, respectively at 5.7 K and 296 K. 

To see these spectra in more detail, we fit the data by 
the sum of two Lorentzian curves La and Lb, 

LM = {^-^Y + ^vl (l = A ' B) ' (1) 

whereby the fitting parameters are the resonant frequen- 
cies, wia and wis, the lincwidths Tia and Tib, and the 
oscillator strengths Sia and Sib- As shown in Fig. 2 
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Fig. 2. (color online) (a) Temperature dependence of the peak 
frequency loia (filled circles) and o^b (open squares). The theo- 
retical curve obtained from the anharmonic potential model 17, 18 
is also plotted as a solid line, (b) Temperature dependence of the 
lincwidth Tia (filled circles) and Tib (open squares), (c) Tem- 
perature dependence of the oscillator strength 5"ia (filled circles) 
and Sib (open squares). The total oscillator strength Sia + Sib 
(crosses) are also plotted. 

(a), o>ia and wib decrease, respectively, by 10% and 
7% from room temperature down to 5.7 K, the latter 
of which closely reproduces our previous results. 10 The 
temperature dependence is well explained by the model 
calculations for anharmonic vibrations. 17 ' 18 On the other 
hand, for the linewidth shown in Fig. 2 (b), Tia de- 
creases monotonically from 0.6 THz at 296 K to 0.4 THz 
at 5.7 K, whereas Tib of 0.06 THz at 296 K (consistent 
with INS measurements 7 ) increases almost by a factor of 
2, to 0.11 THz at 5.7 K. Such an anomalous broadening 
can also be discussed from the difference between fre- 
quencies where o^fw) has a maximum and a minimum 
around 1.2 THz, which is gradually enhanced toward low 
temperature, as seen in Fig 1 (b). It is worth noting here 
that the ratio Sib/ Sia increases from 0.8 at 296 K to 
1.9 at 5.7K as seen in Fig. 2 (c), whereas the total con- 
tribution Sia + >Sib (marked by crosses) is conserved 
within this temperature range. Such a shift of optical 
weights indicates that via and vib are not independent 
vibrational modes, since those do not satisfy the opti- 
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cal sum rule individually. From these readdressed data, 
we conclude that the main sharp peak of the v\ mode is 
superimposed on a broad spectra weighted in the lower 
frequency regime and becomes broadened toward lower 
temperature. 

4. Discussions 

4-1 Model for hybridized phonon system 

Here we propose an impurity-scattering model for 
the hybridized phonon system to interpret the observed 
anomalies in ai(uj). As already introduced, the rattling 
and acoustic phonons are composed of an array of elec- 
tronegative cages and electropositive guest ions. We con- 
sider a simplified model where the cages in the unit cell 
are represented by a single cage at each lattice position, 
as schematically illustrated in Fig. 3 (a). We denote the 
lattice position by X n o, which is the equilibrium position 
of the n-th cage. The mass of the cage and the guest ion 
are Mx and m x , respectively. Then the field variables 
in the phonon system are given by the coordinate of the 
cage, X n Q + X n (the canonical momentum Px n ) and the 
coordinate of the guest ion, X n Q + x n (the canonical mo- 
mentum p xn ). 

Applying an elastic theory to the cage system and a 
local potential model to the guest ion, we consider the 
following Hamiltonian, 



p2 

M X 



J2 ( v <*« 



Kn 



(V • x n 



E 



p 2 

V V(r„) 

TTlrr. 



(2) 



where represents the gradient along the i axis, and 
— X n . By use of a center of mass coordinate, 
R n = (MxX n + m xX n )/(Mx + rn x ) and the relative 
coordinate f n , the Hamiltonian (2) represents both the 
acoustic phonon system with mass M = Mx + m x and 
spring constants K\ and K%, plus the rattling phonon 
system with reduced mass m = m x Mx/(Mx + m x ) in 
the potential V{r). The interaction between R n and r n 
is mediated by 



1 -Tlx 

1 Vx 



(3) 



where r/x = rn x /M and rj x = Mx/M. A strong mixing 
between the rattling and acoustic phonons is expected 
and shown as solid lines in Fig. 3 (b), which is consistent 
with Rcf.. 7 

The local potential V(r) is assumed as 



V ^= E (\ kr 



-Xr 4 



(4) 



for simplicity, as discussed in the one-dimensional an- 
harmonic potential (1D-AP) model. 17,18 We neglect the 
directional difference in the real cage, by assuming the 
averaged effect in the unit cell and concentrating on the 
lowest-lying rattling oscillation. We also assume a quasi- 
on-center rattling as already mentioned. To note, the 
linewidth at 296 K is much smaller than that predicted 
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Fig. 3. (color online) (a) Schematic representation of a rattling 
and acoustic phonon system in the cagelike materials, (b) Zero- 
th order dispersion relation of the phonons propagating along the 
(111) direction, calculated from Eq. (2). The rattling and acoustic 
branches without any hybridization are also plotted in dotted and 
broken lines, respectively . 



from fD-AP. 10,18 The rattling excitation may have a 
more propagating nature through the interaction with 
the acoustic phonon, and the local potential may con- 
tribute rather as a mean field, the energy-spread of the 
level transitions being washed away. The mean field ap- 
proximation in Ref., 

mb&p * k + 3\(rZ), (5) 

and the local potential approximation in Ref. 18 give the 
same result of frequency softening in the present parame- 
ter range. The theoretical curve in Fig. 2 (a) is the model 
result with parameters (1 / '2k) y/kjm = 1.115 THz and 
hX/\/mk 3 = 1.17 x 1CP 2 . The higher order contributions 
from A are negligible, so the Hamiltonian (2) gives no 
information on the linewidth. 

The observed temperature dependence of the spectral 
linewidth motivates the identification of an origin for the 
anomaly. It does not show any increase with tempera- 
ture, which would be expected from certain dissipative 
mechanisms into a continuum of multi-bosonic particle 
states. As previously mentioned, the charge carriers can 
have little effect in phonon scattering. The remaining 
possibility could reside in elastic scatterings from impu- 
rities or defects. Let us consider local potentials of rat- 
tling and acoustic phonons from impurities situated at 

Zg + XlQ] 



H] 



= E 



(6) 



where U c and U g are the force constants of impurity po- 
tentials for the cages and the guest ions, respectively. 

Now the phonon retarded Green function, G(u>,q) is 
obtained from 



? - oj 2 (q) - U rr (uj) 
rixuj 2 (q) - U Rr (iL>) 

i] X uJ 2 (q) - n r _R(cj) 



(7) 



For simplicity, we assume the separation of the longitudi 
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nal and transverse components, indicated by the suffixes 
L and T respectively, and a simple cubic lattice giving 



1 



J T,L 



1 



f — -(cosg^a + cos q y a + cos q z a) 



, (8) 

where u T = A /l2A'i/M, lu l = V / 12(K 1 + K 2 )/M and 
a is the lattice constant. The latter assumption about 
i^>t,l(q) is exact for the type-I clathrate structure. The 
self-energies obtained from the interaction (6) n(w)'s are 
renormalized as 11(0) = 0, consistent with translational 
invariancc, and are expressed as 



U Rr (oj) 



U rR (oj) 

ILrr(uj) 



1 



U 2 G xx (id) 
U c U g G xX {u) 



U c U g G Xx (cd) 
U 2 G xx (lu) 



V: 

1 ~Vx 

1 Vx 



(9) 



Here n,j is the density of impurities, and the equiposi- 
tional Green functions for X n , G xx (uj), and x n , G xx (ui), 
etc, are obtained as 



Gxx(oj) 
G xX {uj) 



G Xx {uj) 
G xx (u>) 



(2tt) 2 



d 3 q 



-Vx 

Vx 



G(uJ,q) 



1 

-Vx 



1 



(10) 

where Qb is in the first Brillouin zone. By iterating 
the equations for the phonon propagators (7), the mean 
field (5) and the self energy (9), we can obtain the 
phonon Green functions and the optical conductivity self- 
consistently. Even when there are several cages in the 
unit cell, Eq. (7) is not much modified, if restricted to 
the lowest rattling vibration and the acoustic phonon. 
Details of the formulation will be presented elsewhere. 

4-2 Calculated optical conductivity and discussions 

According to linear-response theory, the optical con- 
ductivity (j(w) is related to the Green function at the T 
point for r n as 



<j[uj) = iui G rr (U!,q = 0), 



(11) 



because only the relative motion between the cage and 
the guest ion is responsible for the induced electric field. 
Here, N is the density of guest modes, q is the effective 
charge, and G rr (uj, q) is an r-r component of G{u>, q). 

For the numerical calculation of <7i(w) for the v\ mode, 
we fix the parameters as follows. By using the sound 
velocities from the INS measurements, 7 the frequency 
of the transverse and longitudinal acoustic modes at 
the zone boundary are estimated as ujt/2tt = 1.35 THz 
and ujl/2tt = 1.85 THz, respectively. For the rattling 
phonons, we adopt uo/2tt = 1.0 THz, and A = 1.4 x 10~ 2 , 
where ido = ^Jk/m and A = h\/m 2 (dQ are the res- 
onant frequency and the dimensionlcss anharmonicity 
parameter, respectively. Taking into account the stoi- 
chiometry of type-I BGG, the mass ratio is assumed as 
m x /M x = 1/3. The coupling constants are niU 2 /K 2 = 
riiUc/K 2 = 1.73 x 10 -3 . With these parameters the 
calculations are compared with the observed spectra as 
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Fig. 4. Spectral shapes of a\(ui) from v\ at 5.7K and 296K. The 
open symbols and solid lines indicate the data and calculations, 
respectively. The calculated self-energy — (l/7r)Imn rr at 5.7K and 
296 K arc also shown in the inset. Details arc described in text. 



shown in Fig. 4. The experimental data, particularly the 
lincwidth broadening toward low temperature and the 
lincshapes, can be well described by our calculations. 

In order to understand the lincwidth broadening at 
low temperature, the imaginary part Imn rr /7r of the 
self-energy is shown in the inset of Fig. 4. As is obvious 
from Eq. (9), the self-energy reflects the density-of-states 
(DOS) of rattling and acoustic phonons. Therefore, the 
rattling phonon spectra have a sideband structure reflect- 
ing the phonon DOS. Since the DOS is determined self- 
consistcntly, this process gives a strong self-interaction. 
One can see from the inset of Fig. 4, a large modification 
of the self-energy contribution specially around wmf- In 
the present model the temperature dependence comes 
only through wmf, which shows softening toward low 
temperature due to the anharmonicity effect. Strongly 
hybridized rattling and acoustic phonons are responsible 
for a large contribution to the DOS at the zone bound- 
ary. Since the frequency of the rattling phonon is low 
enough to hybridize with the acoustic phonon as the INS 
measurement have revealed, 7 a large DOS renormaliza- 
tion occurs around the frequency of the rattling phonon. 
When the temperature is lowered, wmf is softened, which 
bends down the hybridized acoustic branch and increases 
the DOS at the zone boundary. This gives a large dissi- 
pation in the rattling phonon, filling the would-be hy- 
bridization gap and resulting in wider linewidth at low 
temperature. Both the softening due to anharmonicity 
and the strong hybridization between rattling and acous- 
tic phonons are responsible for the lincwidth broadening 
toward low temperature. 

The above discussions indicate that not only the hy- 
bridization of the acoustic and rattling phonons but also 
impurity scattering in the hybridized phonon system 
is significant for the dynamical properties of clathrate 
compounds. As is obvious from Eq. (9), the impurities 
can also dissipate the heat-conducting acoustic phonons 
around wmf strongly through the large DOS. Here, the 
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origins of such impurity scattering will require further in- 
vestigation but one can consider several candidates, such 
as vacancies and/or substitutions of Ge by Ga in the 
framework. The contribution of impurities should also 
be taken into account for the heat transport behavior in 
this system. 

Although we have not discussed any multiphonon scat- 
tering including the Umklapp process, their importance 
should not be neglected for the suppressed thermal con- 
ductivity in the hybridized phonon system. When the 
two phonons q and q' scatter and the momentum of the 
resulting phonons q" is outside the Brillouin zone, it is 
an Umklapp process. Since the Umklapp scattering can 
take place only when phonons have propagation wavevec- 
tor q' and q" around the zone boundary, such a process 
may not significantly affect the rattler modes at q ~ 0. 
Further investigations on the phonon dispersion relation 
and its temperature dependence, based on INS and/or 
Brillouin scattering measurements, will be required. 

5. Conclusion 

We have performed a detailed experiment of the op- 
tical conductivity measurements for type-I BGG having 
p-type carriers by THz-TDS with improved frequency 
resolution and have further clarified the dynamical prop- 
erties of the infrared-active rattling phonons. Ba(2) 2+ 
modes in the oversized cage show softening and anoma- 
lous broadening toward low temperature. They also have 
broad sideband structures at the lower frequency part of 
the sharp main peak. Our discussion, based on an impu- 
rity scattering model for rattling modes that are strongly 
hybridized with cage acoustic phonons, offers and inter- 
pretation whereby the anomalous features of the rattling 
phonon spectra arise from both the softening due to an- 
harmonicity and the strong hybridization between rat- 
tling and acoustic phonons. 

We also point out that the present experiment has 
disclosed that the anharmonic potential for the rattling 
phonons works as a mean field at room temperature due 
to the strong hybridization with the acoustic phonon, 
and the perturbational treatment is justified. This result 
raises an interesting problem in the off-center potential, 
for example in BGS, since a crossover from a mean field 
regime to a degenerate level transition is expected with 
decreasing temperature. In this sense, interplays between 
carrier density and potential type remains as future prob- 
lems in order to clarify the mechanism operating in cage 
systems. 
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Fig. A-l. (color online) (a, b) Time evolution of THz- wave elec- 
tric fields at low temperature obtained from three different sets 
of measurements (a) in present work and (b) in previous work; 10 
vacuum (broken lines, BG), the substrate only (dot-dashed lines, 
Sub.), and the type-I BasGaieGeso (solid lines, BGG) sample 
glued onto the substrate (solid lines). The second pulses recip- 
rocated in a substrate are indicated by arrows, (c) The real part, 
crijW), of complex conductivity at low temperature in the present 
work (filled symbols) compared with that in the previous work 10 
(open symbols). 



Appendix A: Frequency resolution 

The time evolution of electric fields E{t) at low tem- 
perature measured in the present and previous works, 10 
obtained from three different sets of measurements, vac- 
uum (BG), the substrate only (Sub.), and the sample 
glued onto the substrate (BGG), are shown in Fig. A-l 
(a) and (b), respectively. As indicated by arrows, the 
second pulses reciprocated in a substrate appear only in 
E(t) of the substrate and sample, which delay, respec- 
tively, about 40 ps and 10 ps behind the first pulse. To 
obtain the Fourier-transformed spectra, E(t) after the 
second pulses must be zero-filled by multiplying the fol- 
lowing rectangular window function 



9 T (r) 



1 (0 < r < T) 
{otherwise), 



(A-l) 



where T is the effective length of time delay between 
the first and second pulses. Since the resolution func- 
tions, which are Fourier components of Qt(t), are con- 
voluted into the Fourier-transformed spectra, the actual 
frequency resolution is limited by the thickness of a sub- 
strate. Thus, by using a thicker substrate, the frequency 
resolution can be improved. Figure A-l (c) shows the real 
part <ji(uj) of the optical conductivity obtained in the 
present work compared with that previously reported. 10 
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The linewidth of each observed peak in the present work 
is narrower than that in previous one. 

Appendix B: Analysis for interference effect 

The multiply reflected lights in a sample can interfere 
with each other, which often give rise to the spurious pe- 
riodic pattern in the optical conductivity spectra. Thus, 
we obtain the optical conductivity spectra without an in- 
terference effect 20 between multiply reflected lights in the 
present highly dielectric material (Sec also Ref. 10,11 ) by 
the following procedures. First, the complex refractive 
index is numerically determined by using the equation 
for the transmission coefficient t(ui) taking into account 
multiple reflections in the sample and adhesive; 



t(w) 



t , 



exp{i[(f>f + <j) m - uj{df + d g )/c }} 



x {1 - r fm r fv exp(2%) - r ms r gf exp(2i0 m ) 
- r ms rf v exp(2i</>/ + 2i0 m )}~ 1 , 

(B-l) 

where n, is the complex refractive index, Cq the velocity 
of light, (f 



fli+flj 



and 



hiudi/co the phascshift, ti 

nj = arc the complex transmission and reflec- 

tion coefficients, respectively, and thickness di of each 
medium are taken as df ~ d m -C d s . The subscripts v, 
s, m, and / represent vacuum, substrate, adhesive, and 
sample, respectively. Here to note, the second term in 
the denominator of Eq. B-l describes the interference in 
the sample, the third one in the adhesive, and the forth 
one among the sample and adhesive, respectively. For 
numerical analyses of n/, we assumed h m as a constant 
value of n m — 1.66, which was obtained from another 
THz measurement. Then, the complex conductivity a is 
obtained from 

Airier 



1 



(B-2) 



Next, we estimate the remaining interference effect due 
to the very slight misalignment from the ideal configura- 
tion as follows. This effect, which was observed samplc- 
dependently as well in Fig. A-l (c), as the small dip at 
0.8 THz for example, remain weakly in a even though 
the multiple reflections in the sample are taken into ac- 
count by using Eq. B-l. As one can see from Eq. B-l, 
an interference effect appears when the extinction coef- 
ficient Kf caused by phonons or electrons is small. From 
the denominator of Eq. B-l, the remaining interference 
pattern Ak/(w) in Kf of the sample is approximately 
proportional to 

Ak/(w) oc ln[l — i?exp(— 2K/wd//co) cos(2n/w<i//co)], 

(B-3) 

when the refractive index rif of the sample is much 
larger than that n m of the adhesive and k/ <C 1 (where 



R = \hf - l\ 2 /\h f + 1| 2 is the reflectivity). Then, the 
remaining interference pattern Aeri(u;) in <j\ is approx- 
imated to rif Akuj /2tt from Eq B-2, when Ak <C 1 and 
Kf <C 1. Finally, Aai(u) is estimated from 

Acti(o;) ~A (rifU)/2ir) 

x ln[l — Rexp(—2KfUid/c) cos(2nw d/c)]. 

(B-4) 

Taking into account the temperature- and frequency- 
dependent hf, these spurious effects are most consis- 
tently subtracted from the data by adjusting the am- 
plitude parameter A a of about 3. 
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